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We investigate, for two water models displaying a liquid-liquid critical point, the relation between 
changes in dynamic and thermodynamic anomalies arising from the presence of the liquid-liquid 
critical point. We find a correlation between the dynamic fragility transition and the locus of specific 
heat maxima Cp ax ("Widom line") emanating from the critical point. Our findings are consistent 
with a possible relation between the previously hypothesized liquid-liquid phase transition and the 
transition in the dynamics recently observed in neutron scattering experiments on confined water. 
More generally, we argue that this connection between Cp ax and dynamic crossover is not limited 
to the case of water, a hydrogen bond network forming liquid, but is a more general feature of 
crossing the Widom line. Specifically, we also study the Jagla potential, a spherically-symmetric 
two-scale potential known to possess a liquid-liquid critical point, in which the competition between 
two liquid structures is generated by repulsive and attractive ramp interactions. 



I. INTRODUCTION 

By definition, in a first order phase transition, thermo- 
dynamic state functions such as density p and enthalpy 
H change discontinuously as we cool the system along a 
path crossing the equilibrium coexistence line [Fig.^a), 
path 0] . However in a real experiment, this discontinuous 
change may not occur at the coexistence line since a sub- 
stance can remain in a supercooled metastable phase un- 
til a limit of stability (a spinodal) is reached l] [Fig.^b), 
path 0\. 

If the system is cooled isobarically along a path above 
the critical pressure P c [Fig.^b), path a], the state func- 
tions continuously change from the values characteristic 
of a high temperature phase (gas) to those characteris- 
tic of a low temperature phase (liquid). The thermo- 
dynamic response functions which are the derivatives of 
the state functions with respect to temperature [e.g., iso- 
baric heat capacity Cp = (dH/dT)p] have maxima at 
temperatures denoted T max (P). Remarkably these max- 
ima are still prominent far above the critical pressure 
0j 0. S > an d the values of the response functions at 
Tmax(-P) (e.g., Cp ax ) diverge as the critical point is ap- 
proached. The lines of the maxima for different response 
functions asymptotically approach one another as the 
critical point is approached, since all response functions 
become expressible in terms of the correlation length. 
This asymptotic line is sometimes called the Widom line, 
and is often regarded as an extension of the coexistence 
line into the "one-phase region." 

If the system is cooled at constant pressure Pq, and Pq 
is not too far from the critical pressure P c , then there are 
two classes of behavior possible, (i) If Po > P c (path a), 
then experimentally-measured quantities will change dra- 
matically but continuously in the vicinity of the Widom 



line (with huge fluctuations as measured by, e.g., Cp). 
(ii) If Po < P c (path 0), experimentally-measured quan- 
tities will change discontinuously if the coexistence line is 
actually seen. However the coexistence line can be diffi- 
cult to detect in a pure system due to metastability, and 
changes will occur only when the spinodal is approached 
where the gas phase is no longer stable. The changes 
in behavior may include not only static quantities like 
response functions || but also dynamic quantities like 
diffusivity. 

In the case of water — the most important solvent for 
biological function 0,0 — a significant change in dynami- 
cal properties has been suggested to take place in deeply 
supercooled states JM 0, El EH • Unlike other network 
forming materials [l2| |. water behaves as a fragile liq- 
uid in the experimentally accessible window 0, ll3L Il4[ . 
Based on analogies with other network forming liquids 
and with the thermodynamic properties of the amor- 
phous forms of water, it has been suggested that, at 
ambient pressure, liquid water should show a crossover 
between fra gile behavior at high T to strong behavior at 
low T [3, [la, UjJUjl LUiUil m the deep supercooled region 
of the phase diagram below the homogeneous nucleation 
line. This region may contain the hypothesized liquid- 
liquid critical point |2fJ, the terminal point of a line of 
first order liquid-liquid phase transitions. According to 
one current hypothesis, the liquid-liquid critical point 
is the thermodynamic source of all water's anomalies 
[2fll2lll23 . l2^| . This region has been called the "no-man's 
land" because to date no experiments have been able 
to make direct measurements on the bulk liquid phase 
|2l|. Recently the fragility transition in confined water 
was studied experimentally |2J, l25l |26( since nucleation 
can be avoided in confined geometries. Also, a dynamic 
crossover has been associated with the liquid-liquid phase 
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FIG. 1: Schematic phase diagram for the systems discussed in 
this paper, (a) The critical region associated with a liquid-gas 
critical point. Shown are the two features displaying math- 
ematical singularities, the critical point (closed circles) and 
the liquid-gas coexistence line (bold dashed curve) . (b) Same 
as (a) with the addition of the gas-liquid spinodal and the 
Widom line. Along the Widom line, thermodynamic response 
functions have extrema in their T dependence. Path a de- 
notes a path along which the Widom line is crossed, while 
path p denotes a path crossing the coexistence line, (c) A 
hypothetical phase diagram for water of possible relevance 
to the recent neutron scattering experiments by Chen et al. 
[25H2H on confined water. The negatively sloped liquid-liquid 
coexistence line generates a Widom line which extends be- 
yond the critical point, suggesting that water may exhibit a 
fragile-to-strong transition for P < P c (path a), while no dy- 
namic changes will occur above the critical point (path /3). 
(d) A sketch of the P — T phase diagram for the two-scale 
Jagla model. Upon cooling at constant pressure above the 
critical point (path a), the liquid changes, as the path crosses 
the Widom line, from a low density state (characterized by a 
non-glassy Arrhenius dynamics) to a high density state (char- 
acterized by non- Arrhenius dynamics) as the path crosses the 
Widom line. Upon cooling at constant pressure below the 
critical point (path 0), the liquid remains in the LDL phase 
as long as path /3 does not cross the LDL spinodal line. Thus 
one does not expect any dramatic change in the dynamic be- 
havior along the path f3. 



transition in silicon and silica [27], |28j . In this work, we 
offer an interpretation of the dynamic crossover (called 
a fragility transition or fragile-strong transition by many 
authors) in water as arising from crossing the Widom 
line emanating from the hypothesized liquid-liquid crit- 
ical point [22j [Fig. [Jc), path a]. Our thermodynamic 
and structural interpretation of the dynamic crossover 
may not hold for liquids for which the fragile-strong dy- 
namic crossover can be caused by other mechanisms, as 




FIG. 2: The "two-scale" Jagla ramp potential with attractive 
and repulsive ramps. Here Ur = 3.5Uo, Ua = —Uo, a, is the 
hard core diameter, b — 1.72a is the soft core diameter, and 
c = 3a is the long distance cutoff. In the simulation, we use 
a as the unit of length, and Uo as the unit of energy. 



discussed in [21 



II. METHODS 

Using molecular dynamic (MD) simulations, we study 
three models, each of which has a liquid-liquid critical 
point. Two of the models, (the TIP5P [U and the ST2 
|31|) treat the water molecule as a multiple-site rigid 
body, interacting via electrostatic site-site interactions 
complemented by a Lennard- Jones potential. The third 
model is the spherical "two-scale" Jagla potential with 
attractive and repulsive ramps [Fig. [2] which has been 
studied in the context of liquid-liquid phase transitions 
and liquid anomalies [l6l l32|. For all three models, we 
evaluate the loci of maxima of the relevant response func- 
tions, compressibility and specific heat, which coincide 
close to the critical point and give rise to the Widom line. 
We provide evidence that, for all three potentials, a dy- 
namic crossover occurs when the Widom line is crossed. 

Our results for the TIP5P potential are based on MD 
simulations of a system of N = 512 molecules, carried 
out both in the NPT and NVT ensembles using the tech- 
niques described in [33]. For ST2 simulations N = 1728 
molecules are used and all the simulations are carried 
out in NVT ensemble. For the Jag la p otential, dis- 
crete molecular dynamics simulation [32j | implemented 
for N = 1728 particles interacting with step potentials 
[3l is used in both NVT and NVE ensembles. 



III. RESULTS 

Fig. EI a ) shows for TIP5P the relevant portion of the 
P — T phase diagram. A liquid-liquid critical point is 
observed [33l l35l|. from which the Widom line develops. 
The coexistence curve is negatively sloped, so the Clapey- 
ron equation implies that the high-temperature phase 
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FIG. 3: Results for the TIP5P potential, (a) Relevant part of the phase diagram, showing the liquid-liquid critical point C 
at P c » 320 MPa and T c ~ 217 K, the line of isobaric specific heat maxima Cp ax and the line of isothermal compressibility 
maxima 7('™ ax . (b) Arrhenius plot of the diffusion constant D as a function of 1000/T along different isobars. The filled circles 
indicate the temperatures at which the Cp ax line is crossed, (c) Arrhenius plot of D as a function of 1000/T for P = 100 MPa 
(path a). At high temperatures, D can be fit by D ~ {T — Tmct) 1 (dashed line, also shown in the inset), where Tmct ~ 231 K 
and 7 ~ 1.9. At low temperatures the dynamic behavior changes to that of a liquid where D is Arrhenius (solid line), (d) 
Log-log plot of D as a function of T — Tmct for P = 400 MPa (path 0). The behavior of D remains non- Arrhenius for the 
entire temperature range and is consistent with D ~ (T — Tmct) 1 , with Tmct ~ 201 K and 7 ~ 2.5. Note that the power law 
fits for 7 and Tmct are subject to error due to the relatively small ranges of D and T — Tmct- 



is a high-density liquid (HDL) and the low-temperature 
phase is a low-density liquid (LDL). Fig. EJb) shows the 
T dependence of the diffusion coefficient D, evaluated 
from the long time limit of the mean square displace- 
ment along isobars. The isobars crossing the Widom line 
[Fig. Etc), path a] show a clear crossover from a non- 
Arrhcnius behavior at high T [which can be well fit by a 
power law function D ~ (T — Tmct) 7 ], consistent with 
the mode coupling theory predictions H3)> to an Ar- 
rhenius behavior at low T [which can be described by 
D ~ exp(— E a /T)\. The crossover between these two 
functional forms takes place when crossing the Widom 
line. 

For paths (3 [Fig.EJd)], crystallization occurs in TIP5P 
|33| , so the hypothesis that there is no fragility transition 



cannot be checked at low temperature. Hence we con- 
sider a related potential, ST2, for which crystallization 
is absent within the time scale of the simulation. Sim- 
ulation details are described in [3^1 - This potential also 
displays a liquid-liquid critical point pfl l3?l | , as seen in 
the phase diagram of Fig. 0t a )- The analog of Fig. Etb) 
is shown in Fig. Etb). We confirm that along paths a 
a fragility transition takes place [Fig. Etc)]. Moreover, 
along paths (3 the T dependence of D does not show any 
sign of crossover to Arrhenius behavior and the fragile 
behavior is retained down to the lowest studied temper- 
ature (note that 10 3 /T extends to A.&K^ 1 ). Indeed, for 
paths [3, the entire T dependence can be fit by a power 
law (T- Tmct) 7 [Fig- Etc)]. 

Thus we see that the simulations for both T1P5P and 
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FIG. 4: Analog of Fig. EJfor the ST2 potential, (a) Relevant part of the phase diagram, showing the liquid-liquid critical point 
C at P c « 186 MPa and T c « 246 K, the line of isobaric specific heat maxima Cp ax , the line of isothermal compressibility 
maxima K^"*, and the spinodal lines, (b) Arrhenius plot of the diffusion constant D as a function of 1000/T along different 
isobars. The filled circles indicate the temperatures at which the Cp ax line is crossed, (c) Arrhenius plot of D as a function of 
1000/T for P — MPa (path a). At high temperatures, D can be fit by D ~ (T — Tmct) 1 (dashed line, also shown in the 
inset) where Tmct ~ 268 K and 7 » 1.34. At low temperatures the dynamic behavior changes to that of a liquid where D is 
Arrhenius (solid line), (d) Log- log plot of D as a function of T — Tmct for P = 200 MPa (path 0). The behavior of D remains 
non-Arrhenius for the entire temperature range and is consistent with D ~ (T — Tmct) 1 , with Tmct ~ 217 K and 7 ~ 1.7. 
Note that the power law fits for 7 and T — Tmct are subject to error due to the relatively small ranges of D and T — Tmct- 



ST2 water models support the connection between the 
Widom line and the dynamic fragility transition. It 
is natural to ask which features of the water molecu- 
lar potential are responsible for the properties of wa- 
ter discussed here, especially because water's unusual 
properties are shared by several other liquids whose 
inter-molecular potential has two energy (length) scales 
[27L l28| . We next investigate the two-scale spherically 
symmetric Jagla potential. The Jagla model displays — 
without the need to supercool — a liquid-liquid coexis- 
tence line which, unlike water, has a positive slope, im- 
plying that the Widom line is now crossed along a paths 
with P > P c [Figs.^d) andJSJa)]. There is a crossover 
in the behavior of D(T) when the C™ ax line is crossed 
[Figs. OJb) and l^c)]. At high temperature, D exhibits 
an Arrhenius behavior [Figs.|SJb) and E(c)], while at low 



temperature it follows a non-Arrhenius behavior, consis- 
tent with a power law. Along a f} path (P < P c ), D(T) 
follows the Arrhenius behavior over the entire studied 
temperature range [Figs. E£b) andEJd)]. Thus, the dy- 
namic crossover coincides with the location of the Cp ax 
line, extending the conclusion of the TIP5P and ST2 po- 
tentials to a general two-scale spherically symmetric po- 
tential. 



IV. DISCUSSION AND SUMMARY 

Before concluding, we note that our findings are con- 
sistent with the possibility that the observed dynamic 
crossover along path a is related to the behavior of Cp, 
suggesting that enthalpy or entropy fluctuations may 



5 



0.3 



0.2 



0.1 



(a) Jagla Potential 

max 


^ „ max 

K T 


- Spinodal J* 





0.3 



0.4 



0.5 



10 



10 l 



10" 



10"' 



10"' 



(b) 



A 


P=0.175 




P=0.200 


< 


P=0.225 


• 


P=0.250 


® 


~ max 

°p 




P=0.275 




n max 
°P 


► 


P=0.300 




~ max 










1/T 



10 



10 l 



10" 



10"' 



10 



(C) P>P r (path a) 10 »L 



1.0 




P=0.250 



MCT 



Arrhenius fit 
power law fit 

C ™ 



10 



2.0 



3.0 

1/T 



4.0 



5.0 



10 l 



10 



(d) P<P c (path p) 



1.0 



P=0.225 




2.0 



3.0 



1/T 



FIG. 5: Analog of Figs.[3]and 0|for the two-scale Jagla potential, (a) Phase diagram in the vicinity of the liquid-liquid critical 
point C located at P c ~ 0.24 and T c ~ 0.37, the line of isobaric specific heat maxima C*p ax , the line of isothermal compressibility 
maxima 7f™ ax , and the spinodal lines, (b) The T-dependence of diffusivity along constant pressure paths. Several paths a and 
paths f3 are shown: (i) P = 0.175, 0.200, 0.225 < P c (paths f3 in Fig.^d), along which the system remains in the LDL phase), 
(ii) P = 0.250,0.275,0.30 > P c (paths a in Fig. 0d), along which the system does not remain in the LDL-like state, but the 
dynamic behavior changes from Arrhenius to non- Arrhenius) . (c) D as a function of 1/T for P — 0.250 (path a). At high 
temperatures, the fit is Arrhenius D ~ exp( — 1.59/T) (solid line), while at low temperatures the results are consistent with 
D ~ (T - Tmct) 1 with Tmct ~ 0.27 and 7 ^ 2.7 (dashed line, also shown in the inset), (d) For P = 0.225 (path /?), D is 
Arrhenius for the entire temperature range and can be fit by D ~ exp( — 1.62/T). The unit of D is ayUo/m and the unit of P 
is Uo/a 3 . 



have a strong influence on the dynamic properties. The 
role of Cp is consistent with expectations based on the 
Adam-Gibbs [38| interpretation of the water dynamics 
[Hill and of the fragility transition [Toll27|. 

For both water and the Jagla model, crossing the 
Widom line is associated with a change in the In- 
dependence of the dynamics. In the case of water, 
D(T) changes from non- Arrhenius ("fragile") to Arrhe- 
nius ("strong") behavior, while the structural and ther- 
modynamic properties change from those of HDL to 
those of LDL. For the Jagla potential, due to the positive 
slope of the Widom line, D(T) changes from Arrhenius to 
non-Arrhenius while the structural and thermodynamic 
properties change from those of LDL to those of HDL. 



In summary, our results for water are consistent with 
the experimental observation in confined water of (i) a 
fragility transition for P < P c |26| , and (ii) a peak 
in Cp upon cooling water at atmospheric pressure |4l| . 
Thus our work offers a plausible interpretation of the re- 
sults of Ref. |2fJ| consistent with the existence of a liquid- 
liquid critical point located in the no-man's land. 

We thank C. A. Angell, G. Franzese, J. M. H. Lev- 
elt Sengers, L. Liu, M. Mazza, S. Sastry, F. W. Starr, 
B. Widom, and Z. Yan for helpful discussions and NSF 
grant CHE 0096892 and MIUR-FIRB for support. We 
thank D. Chandler and J. P. Garrahan for helpful crit- 
icisms on the manuscript. We also thank the Boston 
University Computation Center, Yeshiva University, and 



6 



StFX hpcLAB (high performance computing laboratory) for allocation of CPU time. 



[1] Debenedetti, P. G. (1996) Metastable Liquids: Concepts 
and Principles (Princeton U. Press, Princeton). 

[2] Levelt, J. M. H. (1958) Measurements of the compress- 
ibility of argon in the gaseous and liquid phase, Ph.D. 
Thesis (University of Amsterdam, Van Gorkum & Co., 
Assen) . 

[3] Michels, A., Levelt, J. M. H. & Wolkers, G. (1958) Phys- 
ica 24, 769-794. 

[4] Michels, A., Levelt, J. M. H. & de Graaff, W. (1958) 
Physica 24, 659-671. 

[5] Anisimov, M. A., Sengers, J. V. & Levelt Sengers, J. M. 
H. (2004) in Aqueous System at Elevated Temperatures 
and Pressures: Physical Chemistry in Water, Stream and 
Hydrothermal Solutions, eds. Palmer, D. A., Fernandez- 
Prini, R. & Harvey, A. H. (Elsevier, Amsterdam), pp. 29- 
71. 

[6] Robinson, G. W., Zhu, S. -B., Singh, S. & Evans, M. W. 
(1996) Water in Biology, Chemistry, and Physics: Ex- 
perimental Overviews and Computational Methodologies 
(World Scientific, Singapore). 

[7] Hydration Processes in Biology: Theoretical and Exper- 
imental Approaches (1999) ed. Bellissent-Funel, M.-C. 
(IOS Press, Amsterdam). 

[8] Angell, C. A. (1993) J. Phys. Chem. 97, 6339-6341. 

[9] Debenedetti, P. G. (2003) J. Phys.: Condens. Matter 15, 
R1669-R1726. 

[10] Starr, F. W., Angell, C. A. & Stanley, H. E. (2003) Phys- 
ica A 323, 51-66. 

[11] Angell, C. A. (2004) Ann. Rev. Phys. Chem. 55, 559-583. 

[12] Horbach, J. & Kob, W. (1999) Phys. Rev. B 60, 3169- 
3181. 

[13] Prielmeier, F. X., Lang, E. W, Speedy, R. J. & 
Liidemann, H. D. (1987) Phys. Rev. Lett. 59, 1128-1131. 

[14] Lang, E. W. & Liidemann, H. D. (2004) Angew Chem. 
Intl. Ed. Engl. 21, 315-329. 

[15] Ito, K., Moynihan, C. T. & Angell, C. A. (1999) Nature 
398, 492-495. 

[16] Jagla, E. A. (1999) J. Chem. Phys. Ill, 8980-8986. 
[17] Jagla, E. A. (1999) J. Phys.: Condens. Matter 11, 

10251-10258. 
[18] Jagla, E. A. (2001) Phys. Rev. E 63, 061509. 
[19] Tanaka, H. (2003) J. Phys.: Condens. Matter 15, L703- 

L711 . 



[20] Poole, P. H., Sciortino, F., Essmann, U. & Stanley, H. E. 

(1992) Nature (London) 360, 324-328. 
[21] Mishima, O. & Stanley, H. E. (1998) Nature 396, 329- 

335. 

[22] Franzese, G., Malescio, G., Skibinsky, A., Buldyrev, S. 
V. & Stanley, H. E. (2001) Nature 409, 692-695. 

[23] Sciortino, F., La Nave, E., & Tartaglia, P. (2003) Phys. 
Rev. Lett. 91, 155701. 

[24] Bergman, R. & Swenson, J. (2000) Nature 403, 283-286. 

[25] Faraone, A., Liu, L, Mou, C.-Y, Yen, C.-W. & Chen, 
S.-H (2004) J. Chem. Phys. 121, 10843-10846. 

[26] Liu, L., Chen, S.-H., Faraone, A., Yen, C.-W & Mou, 
C.-Y. (2005) Phys. Rev. Lett. 95, 117802. 

[27] Saika-Voivod, I., Poole, P. H. & Sciortino, F. (2001) Na- 
ture 412, 514-517. 

[28] Sastry, S. & Angell, C. A. (2003) Nature Materials 2, 
739-743. 

[29] Garrahan, J. P. & Chandler, D. (2003) Proc. Natl. Acad. 

Set. 100, 9710-9714. 
[30] Jorgensen, W. L. (1982) J. Chem. Phys. 77, 4156-4163. 
[31] Stillinger, F. H. & Rahman, A. (1972) J. Chem. Phys. 

57, 1281-1292. 

[32] Kumar, P., Buldyrev, S. V., Sciortino, F., Zaccarelli, E. 

& Stanley, H. E. (2005) Phys. Rev. E 72, 021501. 
[33] Yamada, M., Mossa, S., Stanley, H. E. & Sciortino, F 

(2002) Phys. Rev. Lett. 88, 195701. 
[34] Buldyrev, S. V. & Stanley, H. E. (2003) Physica A 330, 

124-129. 

[35] Paschek, D. (2005) Phys. Rev. Lett. 94, 217802. 

[36] Gotze, W. (1991) in Liquids, Freezing, and Glass Tran- 
sition, eds. J. P. Hansen, D. Levesque, & J. Zinn- Justin 
(North- Holland, Amsterdam), pp. 287-503. 

[37] Poole, P. H., Saika-Voivod, I. & Sciortino, F. (2005) J. 
Phys.: Condens. Matter (in press). 

[38] Adams, G. & Gibbs, J. H. (1965) J. Chem. Phys. 43, 
139-146. 

[39] Angell, C. A., Finch, E. D., Woolf, L. A. & Bach, P. 

(1976) J. Chem. Phys. 65, 3063-3066. 
[40] Scala, A., Starr, F. W., La Nave, E., Sciortino, F. & 

Stanley, H. E. (2000) Nature 406, 166-169. 
[41] Maruyama, S., Wakabayashi, K. & Oguni, M. (2004) A IP 

Conf. Proc. 708, 675-676. 



